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In numerical simulations of nonabelian plasma instabilities in the hard-loop approximation, 
a turbulent spectrum has been observed that is characterized by a phase-space density of 
particles n(p) ~ p~ v with exponent v ~ 2, which is larger than expected from relativistic 
2 <-> 2 scatterings. Using the approach of Zakharov, L'vov and Falkovich, we analyse possible 
Kolmogorov coefficients for relativistic (m > 4)-particle processes, which give at most v = 5/3 
perturbatively for an energy cascade. We discuss nonperturbative scenarios which lead to 
larger values. As an extreme limit we find the result v = 5 generically in an inherently 
nonperturbative effective field theory situation, which coincides with results obtained by 
Berges et al. in large- TV scalar field theory. If we instead assume that scaling behavior is 
determined by Schwinger-Dyson resummations such that the different scaling of bare and 
dressed vertices matters, we find that intermediate values are possible. We present one simple 
scenario which would single out v = 2. 



I. INTRODUCTION 

Turbulence in nonequilibrium field theory [lj has been studied intensively in scalar field theories 
and more recently in the context of quantum chromodynamics (QCD) @-@]. The nonequi- 
librium dynamics of QCD is of particular interest for the early evolution and thermalization of 
quark-gluon plasma [9] in a heavy ion collision. At least in a weakly coup led quark-gluon plasma, 
anisotropic distributions of particles give rise to plasma instabilities [10| and exponential growth 
of collective fields 11 - . These instabilities provide an infra-red source of energy which can cas- 



cade to higher energies by nonabelian interactions once the corresponding field modes have grown 
to nonperturbatively large amplitudes. For a stationary anisotropic plasma, as considered in the 
hard-loop framework, a quasi-stationary state with scale independent transport to the ultra-violet 
can form. This state is characterized by a linear growth of energy 14- 3], with the formation of a 
power-law spectrum 0, E3] indicative of the appearance of Kolmogorov turbulence [H, 0] . 

In this paper, we employ the approach of Zakharov, L'vov and Falkovich [l|, Q] to find possible 
Kolmogorov exponents in a relativistic gauge theory with an approximately isotropic power law 
spectrum: 

1 < n(p) ~ p~~ v for p* <p < A(t) , (1) 
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which for v > 1 represents a cascade from soft perturbatively unstable modes at scales below 
p* towards the ultra-violet. Distributions of this form with v ~ 2 have been observed in hard- 
loop real-time lattice simulations of an anisotropic purely gluonic plasma 0, G3] • Using classical- 
statistical lattice simulations which do not allow for a separation of scales of hard particles and soft 
collective modes, Ref. [1] has found a significantly lower value at large times, v w 4/3, indicative of 
perturbative Kolmogorov scaling. However, these latter simulations may not yet have probed the 
parametrically separated infrared regime responsible for the different behavior seen in the hard-loop 
simulations. 

Taking into account the renormalizable classically scale-invariant interactions of QCD, we obtain 
the results v = 5/3 and v = 4/3 for energy and particle cascades, respectively, which are familiar 
from scalar cp theory for 2 -R- 2 processes. However it is well known that in the infrared regime 
resummations can play an important role. This is especiall y t rue in gauge theories, where all- 
order resummations lead to effectively non-local interactions [la . In addition to perturbative 
scaling exponents, there may also exist nonperturbative ones which dominate at lower momentum 
scales. We therefore consider the effect of nonperturbative propagator and vertex resummations 
that correspond to a modified scaling law for dressed quantities. 

If we determine this modified law by assuming the same scaling at all loop orders, we obtain 
a Kolmogorov exponent v = 5 for an energy cascade generically in an inherently nonperturbative 
effective field theory situation. This coincides with the nonperturbative infrared fixed point found 
previously by Berges et al. in large- N scalar field theory , but not observed in any numerical 
simulation of nonabelian gauge theory so far. Assuming instead that the scaling behavior of full 
vertex functions is determined by Schwinger-Dyson resummations such that the different scaling of 
bare and dressed vertices matters, we find a range of possible intermediate values, and we present 
a scenario which yields v = 2. 



II. ZAKHAROV TRANSFORMATIONS 

We can formulate the condition for stationary, scale-independent transport in terms of detailed 
balance, by requiring that 'gain' and 'loss' terms balance each other, which means that the rate 
at which energy goes into the cascade is equal to the rate at which it flows out. In Appendix [Al 
we show that this condition is equivalent to the requirement that flux be scale invariant (which is 
the approach used in Ref. [2]). The condition for detailed balance can be written using the closed 
time path (CTP) formalism (see e.g. Ref. [20(). We make the following definitions: 

G ret (x, y) = -iG(x - yo)Mx)>(y)]> = G n (x, y) - G 12 (x, y) , (2) 
G adv (x,y) = iQ(y - x ){[<t>(x),<t>(y)]) = G n (x,y) - G 21 (x,y) , 
G s y m (x,y) = -*({0(x),0(y)}> = G n (x, y) + G 22 (x, y) , 
p(x,y) = mx)A{y)\) = i(G ret (x,y) - G adv (x,y)) = id(x,y) . 
We use analogous expressions for the self energies: 

U ret (x,y) = Un(x,y)+U 12 (x,y) , U adv (x, y) = n u (i, y) + U 21 (x, y) , (3) 
U d (x,y) = U ret (x,y)-U adv (x,y) , H^ m (x, y) = U u (x, y) + U 22 (x, y) . 
Using these definitions, the condition of stationary transport can be written: 
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Z = 9(p ) J d 3 p p{P) ■ L ■ ^) = (4) 

-»■ 0(j>q) J d 3 p (U sym (P)d(P) - U d (P)G sym {P)) = . 

Throughout this paper we use a shorthand notation a function of po and p: f(P) : = f(po,p)- 

All momenta in the cascade are considered soft relative to the ultra-violet end of the cascade. 
Consequently, we can try to solve (jl]) by looking for scaling solutions that correspond to infra-red 
fixed points. In a perturbative regime of QCD, the individual building blocks of a given amplitude 
scale according to: 

G(XP) = \X\- 2 G(P) for G =Re [G ret ] or Re [G adv ] , (5) 
d(\P) = sgn(A) \\\- 2 G(P), 

G sym {\P) = \\\- 2 ~ v G sym {P) , n(XP) = \X\- v n(P) , 
U(XP) = sgn (A) |A| U(P) , V(XP) = V(P) , 

where U represents the 3-point vertex and V the 4-point vertex. 



A. 1 <H> 2 processes 



We begin by using in @ a 1-loop self energy diagram which corresponds to a 1 f> 2 process, 
even though this process is kinematically forbidden for massless particles. We use the method in 
Ref. 21] to perform the sums over the CTP indices and obtain: 



Z~e(p ) J d 3 p J dQdL5(P 2 )5(Q 2 )5(L 2 )5 4 (P + Q + L)U(Q,P,L) U(-L,-P,-Q)F(p ,q ,l ) 



F(p , q , lo) = n(p )n(q )n(l ) 



1 sgn (g ) sgn(/ ) 



n(po) n(q Q ) n(l ) 



(6) 



The integral in ([6]) is identically zero, because the product of delta functions has no support. We 
proceed with the calculation to establish our notation. 

Following Ref. [lj] we perform two different transformations and add the results to the original 
expression. The substitutions are (A2 = Po/q'o, A3 = po/l'o)'- 



Uo -> A 2 po, ^0 -> ^I'o, P -> A 2 q, q -> A 2 p', 1 -> A 2 l'} , 
{q ->• X 3 q' , l Q -4 A 3 p , P "4 A 3 1', q ->■ A 3 q', 1 -> A 3 p'} , 

and they correspond respectively to transforming the set of variables: 

{P, Q, L} A 2 {Q', P', L'} , {P, Q, L} A 3 {L', Q', P'} 



(7) 



(8) 



where we have defined P' = (po>p')- ^he integrand in is symmetric under the interchange 
of any two arguments, so we can combine the results from both substitutions with the original 
expression to obtain a result of the form: 



Z~e(p ) f d 3 p' J dQ'dL'(l + S gn(X 2 ) |A 2 | A + +sgn(A 3 ) |A 3 | A )^(P', Q', L') 



(9) 
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The exponent A = 4 — 2v and is obtained from: 

^3-34-1+2— 4|^|— 6+2— 2f _ ^8m-4-2i/ 



1 4- 2i/ 



(10) 



The exponent of the first term is: 3-3 from the integrals over the spatial momenta {p,q, 1}, 1 
from the Iq integral, 2 from the go integral (since go ~~ * X 2 po = Po/q'o) an d —4 from the overall 
momentum conserving delta function. The exponent of the second term is: —6 from the on-shell 
delta functions, 2 from the vertex functions, and — 2u from the distribution functions. For A = — 1 
the bracket in the first line of (O) becomes: 



1 1 1 

l + \- + \- = - 

X 2 A 3 po 



Po + q'o + l'o 



0, 



(11) 



where we have used that the integrand in ([6]) contains a delta function of the form 5(po + q' + 1' ). 

If Z were not identically zero anyway, due to kinematic constraints, A = — 1 would have given 
a scaling solution of v = 5/2. The function F(po,qo, Iq) in (j6]) is zero for v = 1, and therefore v = 1 
is also a solution of Z = 0. This is not a scaling solution however, because the quantity in square 
brackets in is not zero. The solution v = 1 corresponds to the equilibrium thermal distribution 
function, and the equality F(po, qo,lo, ro) \ u=1 = is just the KMS condition. 

B. 2 <H> 2 processes 

Since the contribution to Z from the 1-loop diagram is identically zero due to kinematic con- 
straints, we turn to the 2-loop self energy. If we include all 2-loop self energy diagrams in ([3]), 
corresponding to all possible 2 •<-)■ 2 processes, we obtain: 



Z~0(pq) J d 3 p J dQ dRdL 5(P 2 ) 5{Q 2 ) 5(L 2 )5(R 2 )5\P + Q - L - R)\M\ 2 F(p ,q ,l ,r ) 



F(po,qo,lo,r ) = n(po)n(q )n(l )n(r ) 



1 . sgn(go) sgn(Z ) sgn(r ) 



+ 



n(p ) n(q ) n(l ) n(r ) 



(12) 



where the square of the matrix element \M\ is a symmetric function of the Mandelstam variables: 
s = (P + Q) 2 , t = (P — L) 2 and u = (P — R) 2 . In order to find the scaling solutions to (jU we 
perform three different substitutions and add the results to the original expression: These three 
shifts are: 

Po 

{q ->• X 2 po, Iq -> X 2 l' , r ->• X 2 r' , p -> A 2 q , q -> A 2 p', 1 ->• A 2 l', r ->• A 2 r'} ; A 2 = — , (13) 

% 

{q ->■ -A 3 g , / ->• A 3 p , r -> -A3r , p -> A 3 1', q ->• -A 3 q', 1 -> A 3 p', r -> -A 3 r'} ; A 3 = , 

{go ->• -A 4 g , l o -> A 4 po, -> -A4J0J P A 3r', q -> -A 4 q', 1 -> A 4 p', r -> -A 4 1'} ; A 4 = ^ . 

r o 

Using the notation P' = (po>p') we can rewrite (|13p in the form: 

{P,Q,L,R}^X 2 {Q',P',L',R'}, (14) 
{P, Q, L, R} + A 3 {L', -Q', P', -R'} , 
{P, Q, L, R} -+ A 4 {P', -Q', P', -L'} . 
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Combining the contributions from all three shifts with the original expression, the integral in (|12p 
becomes 1 : 

Z~9(p ) J d 3 p' J dQ'dL'dR' P(P',Q',L',R')[l + sgn (X 2 ) |A 2 | A - sgn (A 3 ) |A 3 | A - sgn (A 4 ) |A 4 | A ] , 

(15) 

with A = 4 — 3v obtained from: 

A 12| A |-8-3, = | A |4-3, = | A |A_ (16) 

The first exponent is (4 • 3 + 4) — 4 from the Jacobian of the transformation and the 4-dimensional 
delta function, and the second exponent receives contributions: —8 from the on-shell delta functions 
and —3v from the distribution functions. 

Using the definitions of the A variables in (I13p together with momentum conservation P + Q = 
L + R, a scaling solution to (JI|) is obtained: 

5 

-l = A = 4-3z/ ^>v = - 

This result is the same in scalar theories when 2 o 2 processes are taken to give the dominant 
contribution to the collision integral. This indicates a universality corresponding to an insensitivity 
to the particular theory under consideration. 

As in section Hi A\ we have an algebraic solution when v = 1. This solution corresponds to the 
equilibrium KMS condition, and is not a scaling solution, because the quantity in square brackets 
in (|15p is not zero. Note also that there is a scaling solution for A = (y = 4/3) if we make use 
the restriction imposed by the delta functions in (|12p to consider only 2 — >• 2 processes, for which 
two of A2, A3 and A4 are positive. This solution corresponds physically to the scale invariance 
of particle flux (see Appendix [A]) . For a relativistic field theory, where particle number is not 
conserved, the physically relevant solution would seem to be the one that corresponds to the scale 
invariance of energy flux. Remarkably however, the numerical results in classical statistical field 
theory simulations of a nonabelian gauge theory in Ref. [8[ appear to favor v = 4/3. On the 
other hand, simulations of nonabelian plasma instabilities within the hard-loop framework have 
found v ~ 2 in a momentum regime that is parametrically separated from the hard scale of plasma 
constituents 0, E3] and perhaps out of reach of the simulations of Ref. . 



(17) 



C. Zakharov Scaling at Arbitrary Loop Order 

We can generalize the results of the previous two sections and obtain an expression for the scaling 
exponent of an /-loop amplitude. In the next section we show how to use this result to find solutions 
to (jl]) that correspond to infinite resummations of vertex corrections. In the CTP formalism, n- 
point functions have 2 n — 1 independent components, and integral equations for n-point functions 
couple different components. However, the situation becomes vastly simplified in the statistical 
limit where we only need to keep terms that carry the greatest number of distribution functions, 
which is equal to the number of loops. We define: I = number of loops, I = number of internal lines, 
v 3 = number of cubic vertices, t>4=number of quartic vertices, and E = number of external legs. We 
denote the complete set of momentum arguments for the external legs: Pi := {Pi, P2, ■ ■ ■ Pe}, and 



1 We use A for the scaling exponent of Z throughout this paper, without introducing subscripts to distinguish results 
that correspond to different choices for the self-energies in (Q. 
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the momentum arguments for the loop variables: Lj := {Li, L2 ■ ■ ■ £/}. An arbitrary amplitude 
can be written symbolically: 

' (G s y m (P u L,)) 1 (G(P i ,Lj)) I ~ l {U{P h Lj)) v * (V(P u L,)r. 

Defining Lj = XL'- and using ([5]) we obtain: 

A{XPi) ={ S gn{\)) E \\\-~<A{P i ), 1 = 2I + v-l-v 3 -Al. (18) 

The exponent on the sign function is obtained from the fact that an amplitude with E odd will 
always correspond to v 3 odd. Using the topological constraints: 

l = I-(v 3 + v 4 ) + l, E = 4t> 4 + 3v 3 - 21 , (19) 

the exponent 7 can be rewritten in different forms: 

7 = v - 4 + E + y z/ + v A v = E - 4 + Iv . (20) 

We can find the Zakharov scaling solution of (dJ with /-loop self-energies by following the method 
in the previous section, provided we choose a set of diagrams that have the appropriate symmetry 
under permutations of the momentum arguments of the on-shell propagators. Note that in the 
2-loop calculation discussed in section III the complete set of 2-loop diagrams produces the fully 
symmetric matrix element in Eq. (I12p . At higher loop orders we expect that the set of diagrams 
with the correct symmetry is equal to the full set of diagrams at that order. The exponent A is 
obtained from —7 by including the two factors (see Eq. (jlj)): 

(1) A 3+1 = |A| 4 where the 3 is from the integral over 3-momenta, and the 1 is from the integral 
over the zeroth momentum component that is interchanged with the external momentum by the 
Zakharov transformation. 

(2) |A|~( 2+i/ ) from the extra factor proportional to G sym (p). 

Combining these factors the solution for the scaling exponent is obtained from: 

-l = A = 4-(2 + i/)-7. (21) 
Substituting ([20|) into ([2Tj) and using E = 2 we obtain: 

5 5 
1+7 ~ l + ±(v 3 + 2v 4 ) ' 

Note that I = 2 gives v = 5/3 as discussed above. Higher loop orders, which correspond to processes 
involving m = I + 2 > 4 particles, are seen to be subdominant. 

In section III Bl we saw that there is also a scaling solution for A = at the 2-loop level 
corresponding to a particle cascade. With arbitrary loop orders contributing, we do not expect 
A = to be a solution. This would require that only n O n processes contribute, which seems 
unlikely in a relativistic field theory, where particle number is not conserved. 

III. RESUMMATIONS 
A. Resummed Vertices 

In section [TT] we have assumed bare 3- and 4- vertices and standard scaling behavior of propaga- 
tors, except for modified statistical factors. When calculating transport coefficients and production 



A(Pi) = 
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rates in gauge theories, the non- locality of effective interactions requires us to include diagrams 



that correspond to infinite series of interactions with soft particles la, LL9|]. Traditionally, these 
diagrams are packaged into resummed vertex corrections. We expect that the Zakharov scaling 
exponents will also receive contributions from processes involving an unlimited number of particles. 
We now consider the possibility that these processes can be included using resummed propagators 
and vertices. 

First we remind the reader of the structure of a traditional ladder resummation. The basic 
idea is that although each term in the ladder series is formally of higher order in the coupling 
constant, kinematic enhancements change the naive power counting. The dominant contribution 
to each ladder diagram comes from the region of phase space that corresponds to hard rails (to 
maximize phase space) and the soft rungs (to obtain a Bose enhancement). The end result is that 
all diagrams in the series contribute at the same order, which means that all ladder diagrams have 
to be resummed to obtain the leading order contribution. 

The structure of the calculation of a scaling exponent is quite different. One basic difference 
is that all lines carry momenta that correspond to soft particles in the cascade, and thus it is not 
possible to separate two scales which could play the role of hard rails and soft rungs. We note that 
the method that we have described to find scaling solutions requires that all distribution functions 
be approximated by their classical forms, so that the integrand satisfies a scaling law. This means 
that, in any scaling calculation, it is not possible to treat some momenta as parametrically harder 
than others, without losing the symmetry that is needed to do the calculation. 

Another difference is discovered by considering the scaling behaviour of an ra-loop amplitude 
calculated with bare lines and vertices. From equation (I18p we see that in the infra-red where |A| 
is large, we have «4(Pj) large for 7 positive and small for 7 negative. From (|20p we have: 

7 tree = £_ 4) y-loqp = E _ 4 + v > ^2) 

and therefore for v positive, a 1-loop vertex calculated with bare lines and vertices is larger than the 
corresponding tree vertex. This situation is completely different from that encountered in ladder 
resummations, where the resummation involves an infinite series of diagrams which all contribute at 
the same order. In our case, higher loop graphs are larger than lower loop ones, and any traditional 
resummation will fail to converge. This is similar to the well known Linde problem, and is therefore 
expected at very soft momentum scales. 

The key point is that in a scaling calculation, we do not need to actually perform the resumma- 
tion, but only to extract effective scaling exponents that correspond to resummed lines and vertices. 
We define a resummed propagator G(P), 3- vertex U{Pi) and 4-vertex V(Pi) with modified scaling 
laws (compare with equation ([5])): 

G(XP) = \X\- 2 - X2U G{P) for G =Re [G ret ] or Re [G adv ] , (23) 

<f(AP)=sgn(A) |A|- 2 -^G(P), 

G syrn {\P) = \\\- 2 -( 1+x ^ G sym {P) , 

U{XP) = sgn (A) |A| 1+ ^ U{P) , V{XP) = \X\ XAV V(P) . 

In addition, we include effective (n > 4)-vertices with scaling laws: 

V n (XPi) = (Sign (A)) n iXr^-^VniPi) . (24) 

The topological relations obtained by generalizing f)19|) are: 

oo oo 

l = I + l-J2 v k, E = Y,kv k -2I . (25) 

k=3 k=3 
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Note that setting x n = we recover the scaling laws of the bare quantities (Eq. ©), apart from 
the fact that there are no bare n-point-vertices with n > 4. 

Equation (|18p gives the scaling exponent for an arbitrary amplitude calculated with bare lines 
and vertices. The generalization of this result for amplitudes with resummed lines and vertices is: 

7x . = 1(2 + x 2 v) + {v- 4)1 - u 3 (l + x 3 v) - f 4 (0 + Xiv) 

(26) 

and using (f25l) to eliminate I and / we find: 

E v 

lx = v-4+—(2-v- x 2 v) + y (1 + 3x 2 - 2x 3 )zy + u 4 (l + 2x 2 - x A )u + ... (27) 

If an all-order nonperturbative resummation is indeed able to change the scaling behavior of an 
arbitrary amplitude, we expect that the scaling exponent of the inverse propagator (a) in Fig. [1] 
equals that of graph (b), and that the 3-point graph (c) has the same scaling exponent as graph 
(d), and that the 4-point graph (e) has the same scaling exponent as graph (f). 



(-i) 




(a) (b) 




(c) (d) (e) {f) 



FIG. 1: Graphs that are required to have the same scaling exponent. The dots indicate any number of 
additional rungs. 



This leads us to the conjecture that a nonperturbative fixed point is characterized by values of 
{x2,X3,X4, . . .} for which 7 only depends on the number of external legs, and not on ^3 and 
and higher. We comment that this includes not only ladder graphs, as shown in Fig. [1] but 
also graphs which correspond to crossed ladders, and even graphs with nondadder structure. As 
discussed previously, the fact that our definition of resummed lines and vertices includes nondadder 
topologies is consistent with the physics of the cascade, in which all particles are soft. Imposing the 
condition that (p7|) is independent of the number of vertices fixes the scaling of effective vertices 
in terms of x 2 according to: 



X4 = 1 + 2x 2 ■ 



(28) 



We could attempt to find the scaling exponent directly from the generalization of (1211) : 



= 3 + l- A-((2 + x 2 v) + v) - 7a 



(29) 



E=2 



Substituting (|27p and (|28|) directly into (i29j) we find that the dependence on x 2 cancels exactly, and 
we obtain v = 5, at any loop order. The cancellation of x 2 can be seen directly using simple power 
counting, and the fact that the integrand in the second line of (Jlj) has the structure n(p)G(p)TL(p) ~ 
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n(p)G(p)(G~ l (p) — Gq (p)). We expect that the result is independent of the number of loops, since 
the resummation procedure we have used is just the condition that diagrams with an arbitrary 
number of loops contribute at the same order. 
The result (in 3+1 dimensions) is: 

v = 5 for A = —1 (energy cascade) . (30) 

This result coincides with the nonperturbative infrared fixed points found in Rcfs. @, @] in (large- 
iV) scalar field theory. Note also that result in (]30p contains a prediction of the values of the Xj's 
which could be checked numerically. As discussed above, the actual scaling of propagators and 
vertices does not play a role, which suggests that this extremely large result is a universal limit. 
The result depends on the implicit assumption of the existence of an effective theory in the infrared 
where exclusively dressed vertex functions govern the dynamics. 

In situations where the specific microscopic dynamics is relevant, the diagrammatics is instead 
governed by Schwinger-Dyson equations which combine bare and dressed vertices. If we assume 
that j x \e=2 is characterised by the diagrams shown in Fig. [2j we should make the substitution 

Vk -> n k + x k vn k (31) 

in (I26D . where n k and n' k are the numbers of bare and dressed A:-vertices, respectively. 



FIG. 2: The Schwinger-Dyson equation for the gluon self-energy. Dressed vertices are indicated by filled 
vertex dots and propagator lines are assumed to be dressed. 

The second diagram in Fig. [2] corresponds to rt3 = n' 3 = 1 and yields: 

forA = -l. (32) 



2 + 3^2 — X3 

Using the scaling relations (|28p , which express nonperturbative behaviour in the sense that all loop 
orders contribute equally to dressed quantities, we obtain: 

LU forA = -l. (33) 



3(1 + x 2 



The tadpole diagram in the Schwinger-Dyson equation does not have an imaginary part. The third 
diagram in Fig. [2] gives u = 5/(3 + 5x2 — 2x3), and the fourth gives v = 5/(3 + 4x2 — £4)- Using 
(f28l) . both of these results reduce to v = 5/ [2(1 +£2)]. Comparing with (j33j) . these diagrams thus 
give subdominant contributions. 



B. Scaling Exponent from Absorbing Propagator Scaling into Resummed Vertices 

As a possible scenario for obtaining a more definite prediction from simple assumptions of non- 
perturbativeness, we add one further conjecture which leads to a different procedure for calculating 
7. This procedure singles out the exponent v = 2 which has been observed in numerical hard- loop 
simulations 0, G3] ■ 
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In (|27p we have used the topological relations ()25[) to write 7^ without reference to the number 
of loops I or internal propagators /. We now assume the effect of a nonperturbative scaling 
of resummed propagators to be equivalent to a nontrivial scaling law of vertices such that it 
corresponds to effectively replacing the topological constraints (|25j) by 

I = ^(2 - E + v 3 ) ->• i[2 - E + v 3 + x 3 vn' 3 ) , 

I = l(-E + 3v 3 ) -+l[-E + 3(v 3 + vx 3 n' 3 )] , (34) 

and that we can then set X2 = 0. We have anticipated the observation (discussed at the end of 
section fill Ap that higher n-point vertices produce subdominant contributions that can be dropped. 
This procedure yields: 

7efr = v - 4 + E[l - -v] + -vv 3 + -vx 3 {{2 + v)n' 3 - 2n 3 ] . (35) 

Having set X2 = 0, (|28p implies x 3 = g, X4 = 1, which means that resummed vertices are 
smaller than bare vertices in the infrared. In ([35]) the terms that contain n 3 and 723 contribute 
with a positive and a negative sign, respectively, which indicates that effectively large resummed 
lines compensate the contribution from small resummed vertices. 

For the second diagram in Fig. El using (|29j) and (l35j) . with X2 = 0, n 3 = n' 3 = 1 and A = — 1 
we obtain: 



7eff| i?= 2 = 1 v = 2, (36) 

in agreement with the result obtained by the numerical simulations in Ref . 0, E3] ■ 

As a consistency check, we can reconstruct the effective resummed propagator implied by the 
substitutions (134p by writing generally: 

7eff = {v~ 4)/ + 7G^o + IgI + lu{v 3 - n 3 ) + jun 3 , (37) 

where Iq and / are the number of bare and resummed lines, respectively, and the parameters 
7G = 2, 7[/ = — 1 and 777 = —2 are the scaling exponents for the bare propagator, bare vertex, and 
resummed vertex, respectively. For the dominant second diagram of Fig. [2] we have: 

1 = 7effU= 2 = ((" - 4) x 1) + ( 76 x 2) + (-1 x 1) + (-2 x 1) ld = 3 . (38) 

This is consistent with the Ward identity K ■ U(-P - K,K,P) = G~ l {P + K) - which 
implies — 1 + 7^ = — 1 — 2 = 7^-1 = —3. 



IV. CONCLUSIONS 

Assuming high occupation numbers such that particle distributions in turbulent cascades can 
be approximated by their classical limit, we have considered perturbative and nonperturbative sce- 
narios for scaling solutions. Perturbative scaling solutions in nonabelian gauge theories correspond 
to 2 — > 2 processes, with exponent v = 5/3 (y = 4/3) for energy (particle) cascades. 

We have also done a non-perturbative analysis and obtained two results using different senarios. 
These results must correspond to different, parametrically separated infra-red regimes. This is at 
least consistent with recent numerical calculations in which different authors have produced both 
results using methods which may probe different regions of the full momentum phase space. In the 
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case of maximal nonperturbativeness, namely an effective field theory where the scaling exponent 
of vertex functions does not become subleading at higher loop orders but has all loop orders 
contributing equally, we obtained the (model independent) result of v = 5 for energy cascades. This 
coincides with the nonperturbative infrared fixed point obtained for large- iV scalar field theories of 
Refs. [3, .4 which is however too large to match existing numerical simulations of nonequilibrium 
nonabelian gauge theories. 

By considering Schwinger-Dyson equations for full vertex functions, where bare and dressed 
vertices contribute separately with perturbative and nonperturbative scaling exponents, we have 
found that intermediate values for v appear possible. We have proposed a scenario for fixing these 
scaling exponents in a simple way that yields v = 2, consistent with numerical results from hard- 
loop simulations of nonabelian plasma instabilities at nonperturbatively large field amplitudes. 



Appendix A: Scale Invariance of the Flux 

We have calculated the exponent v from the requirement of stationary transport, by finding a 
scaling solution to equation (JH) . In this appendix we show that (j3]) is equivalent to the requirement 
that the energy flux be scale invariant, which is the approach used in Ref. 0]. We calculate the 
energy flux through a sphere of radius A and require that the result be independent of A. This 
flux is given by: 

S~ f A d 3 p-p-I(p), J(p) = ^M. (Al) 



Jo dt 
We define the scaling exponent of the collision integral: 

I(\p):=\~ e I(p). (A2) 

To find the momentum dependence of the collision integral we choose A ~ 1/p and write I(p) = 
p 1(1). Integrating we find: 

S(P) ~ P T \ . T = 4-9. (A3) 

For scale invariant flux we require: 

r = and lim = finite , (A4) 

t— >o r 

which means that the collision integral must have a zero at r=0. 
Using equations (jl]) and (|Aip we have: 

Z = 6(p ) J d 3 p p(P) P I(p), (A5) 

and therefore it appears that if I(p, t) has a zero at some value of v then so does Z . However, 
while this is true for algebraic solutions, like the thermal equilibrium solution v = 1, it does not 
mean that a scaling solution of Z = is a zero of the collision integral. To see this note that if we 
multiply the integrand in (|12p by a function f'(P, Q, L) that is symmetric under the interchange 
of any two variables, we would obtain a different critical exponent. From (|Aip and (|A5p we obtain 
that the scaling solution of Z = is given by: 

-1 = A = -2 + 1 + t. (A6) 
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The -2 is the scaling dimension of the factor p(P), and the 1 is the extra factor introduced by the 
Zakharov transformation (see the discussion under Eq. ([!])). From (lA6j) we have that A = —1 is 
equivalent to t=0, which shows that the conditions of stationary transport and scale invariant flux 
are equivalent. 

We note that we would also obtain a solution for scale invariant flux if we remove the factor 
(p) in equation (jAlh . Physically, this corresponds to the scale invariance of particle flux. From 
equation (|A6j) it appears that r — > r — 1 indicates a scaling solution that corresponds to A = 
corresponding to a particle cascade. However, as discussed in section III C\ we do not find the 
A = scaling solution in a nonperturbative scenario unless we add the assumption that only 
particle number conserving processes n o n processes contribute. 
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